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ABSTRACT 

We derive an analytical solution to the computation of the output of a Lyot coro- 
nagraph for a given complex amplitude on the pupil plane. This solution, which does 
not require any simplifying assumption, relies on an expansion of the entrance complex 
amplitude on a Zernike base. According to this framework, the main contribution of 
the paper is the expression of the response of the coronagraph to a single base function. 
This result is illustrated by a computer simulation which describes the classical effect 
of propagation of a tip-tilt error in a coronagraph. 

Subject headings: instrumentation: adaptive optics - techniques: high angular resolu- 
tion 
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1. Introduction 

The discovery of extrasolar planets is at the origin of a renewed interest in stellar coronagraphy. 
Considering the ambition of the targeted objectives, many authors have pointed out the necessity 
for a very accurate analysis of the system in order to study various undesired effects. For example, 
the specific properties of the light intensity measured by a system based on an extreme adaptive 
optics system and a coronagraph are the result that neither the residuals of the turbulence, nor the 
ideal coronagraphed point-s pread function can be neglected with respect to the faint object (planet). 
Aime and Soummer I ( 2004 ) analyzed the fact that the wavefront amplitudes associated to these two 
contributions will interfere leading t o the so-called "pinned" speckles. Another example is given by 
Llovd and Sivaramakrishnan I (j2005l ) which pointed out that a small misalignment of th e star with 



the c enter of the stop can result in a fake source. A related problem is also present in ([Soummer 



2004 ) wich derives the opti mal apodization for an arb itrary shaped aperture using an algorithm 
proposed independently in (jGuvon and Roddier Il200d ) which relies on iterated simulations of the 
coronagraph response. 

More generally, an intense activity aims to optimize the different coronagraph parameters 
(mask size, apodization shape,...) for a number of projects dedicated to devise high-dynamic range 
i maging on the VLT (S phere), Gemini (GPI) or the Subaru telescope (HiCIAO), see for example 
(jAime and Vakili II2006I ). The input/output relation of a coronagraph is in this case simulated by 
numerical computations based on discrete Fourier transforms. However, such a numerical technique 
suffers from the well-known problems related to the choice of the extent of the sampled surface and 
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the sampling frequency which both define the sampling in the transformed domain. Note that this 
compromise is coupled with the difficulty to evaluate numerically the simulation errors. 

This work focuses on the analytical characterization of the response of a Lyot coronagraph. 
The objective is obviously also to gain deeper insight in the behaviour of the system. This prob- 
lem has already been studied in the literat ure and analytical results were obtai ned under various 
assumptions. In the one-dimensional case, iLlovd and Sivaramakrishnan I (|2005l ) assume that the 
Lyot stop is b and limited and the phase on t he telescope aperture is small. This last hypothesis 



is removed in (jSivaramakrishnan et al. II2005I ) where the computation is carried for a rectangular 



pupil assuming again that the Lyot stop is band-limited. The development presented herein for 
a circular pupil differs from these approaches substituting these simplifying assumptions by an 
expansion of the complex amplitude on an orthogonal basis. 

Section 1 recalls the general formalism of Lyot coronagraphy and justifies the choice of an 
expansion of the complex amplitude on a Zernike base. Section 2 contains the main results of 
the paper; the response of the coronagraph to a Zernike polynomial is computed. The result 
involving an infinite sum, a bound on the truncation error is then derived. Section 3 presents two 
simulations. First the response of the coronagraph to the 6 first Zernike functions is computed. 
Then the formalism derived in this paper is used to illustrate the effect of a tip-tilt error in a 
coronagraph. A short appendix containing the material required for the mathematical derivations 
of section 2 is included at the end of the paper. 



2. Notations and hypothesis 



2.1. Coronagraph formalism 



We follow the notations of Aime et al. ( 2002i ) and Soummer et al. ( 2003 ). The successive 
planes of the coronagraph are denoted by A, B, C and D. A is the entrance aperture, B denotes 
the focal plane with the mask (without loss of generality we assume that the amplitude of the 
mask is 1 — e where e = 1 corresponds to the classical Lyot coronagraph and e = 2 to the Roddier 
coronagraph), C is the image of the aperture with the Lyot stop and D is the image in the focal plane 
after the coronagraph. The aperture transmission function is p{x, y) and the wavefront complex 
amplitude in A is ^'(x, y). In the case of an apodized pupil, we assume that the apodization function 
is included in ^'(x,y). In order to simplify the notations, the mask function in B is defined with 
coordinates proportional to 1/A/ and decomposed as: 



1 — em — — , -— 



where function m(.) equals to 1 inside the coronagraphic mask and outside. 



(1) 



We will make in the sequel the usual approximations of paraxial optics. Moreover we neglect 
the quadratic phase terms associated with the propagation of the waves or assume that the optical 
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layout is properly designed to cancel it (jAime 1 120031 ) . The expression in cartesian coordinates of 
the complex amplitude in the successive planes are: 



(^-aC-x, -y) - e [i!A{-u, -v) * in {u, v)] (x, y)) p{-x, -y) 



1 ( X y . , , 
-.^B ( — ' — I VK-x, -y) 



(2) 
(3) 

(4) 
(5) 

(6) 



where / is the Fourier transform of / and * denotes convolution. Eqs. (j5l6p assume that the Lyot 
stop is the same as the pupil. However for classical "unapodized" Lyot coronagraph the residual 
intensity in plane C is concentrated at the edges of the pupil and a reduction of the Lyot stop size 
is needed in order to improve the rejection. The case of a reduced Lyot stop, which consists in 
convolving Eq. ([6]) by the appropriate function, has not been considered in Eqs. (j5|6p to alleviate 
the notations but will be discussed in section 3. It is important to note that the reduction of the 
Lyot stop can be avoided using a prolate apo dized entrance pup il which will optimally concentrate 



(*A(-a;, -y)m{-x, -y) * p{-x, -y)j jj 



the residual amplitude in C, see for example (jAime et al. II2002I ). 



The coronagraph response being derived herein for a circular pupil, the use of polar coordinates 
will be preferred. Transcription of previous equations to polar coordinates is straightforward. 
Moreover, as long as the aperture transmission function and the stop have a circular symmetry, 
their Fourier transform will verify the same symmetry, as proved by Eq. (lASh with m = 0, i.e. the 
Hankel transform. This leads to the following expression of the complex amplitude in D: 

-1 



-^A(r,^ + vr) + — 
jA/ jXf 



^'yi(r, 9 + 7r)m(r) I * p{r) I (r, 9) 



(7) 



where the convolution of the two functions is still computed with respect to to the cartesian coor- 
dinates (x, y). 



2.2. Choice of a base 

As mentioned in the introduction, the analytical computation of the coronagraph response 
proposed herein relies on the expansion of the complex amplitude in A on an orthogonal basis. 
Eq. d?]) shows that the coronagraph acts linearly on the complex amplitude, consequently the 
problems simplifies to the computation of the response of each basis function. The retained solution 
consists in the expansion of the complex amplitude in A on Zernike polynomials. Basic properties 
of the Zernike polynomials required in the paper are recalled in appendix lAl 
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Adopting the usual ordering of the Zernike circle polynomial ([Mahaian Ill994l ) we can write: 

(8) 



^Air,e) = Yl a^m,n)U::'{r/R,9) 

{m,n) 



(9) 



where R is the radius of the aperture. This expansion is rather unusual, the Zernike polynomials 
being generally used for the expansion of the wavefront. However it is worthy to note that, as 
Eq. ([5]) shows, a coronagraphic system will always introduce amplitude aberration. Hence, even in 
the case of a perfect wave with no aberration in A, an expansion of only the phase in C will not 
be appropriate. Finally, Eq. ^ can also be justified by the fact that it coincides (up to a linear 
transform) with the classical approximation of the complex amplitude in the case of sufficiently 
small phase errors assuming a first order development of the exponential function. 

We will illustrate the expansion ([9]) in the case of tip-tilt error with an apodized pupil: 

■^A{rR, 0) = a(r)n(r)e^^''"°'^('') (10) 

where a{r) denotes the pupil apodization and n(r) = 1 for r G [0, 1) and if r > 1. Computation 
of the projection of ^/^(r, ^) on U^(r/R,9) is straightfo rward using the definition of the Bessel 
functions of integer order (jAbramowitz and Stegun Ill972l ): 

r.27r f-R fl j-2-K 

'^Air,9)U;^{r/R,e)rdrde = R^ / R';^{r) cos{me)a{r)e^^''''''''^^^rdrde (11) 
Jo Jo Jo 

1 

a{r)R'^ir)Jm{Pr)rdr (12) 





The projection of ^'^(r, 6*) on U^^"^{r / R,6) equals 0. 

• In the unapodized case, a(r) = 1, integral in Eq. (fT2]) can be computed using Eq. \A2\ : 



27TR^r C a{r)Rl\r)Jm{Mrdr = 27ri?2j™(-l)^ 
Jo 



/3 



(13) 



The coefficie nt is then obtaine d dividing this quantity by the norm of the Zernike 
polynomials (jBorn and Wolf 1119911 ). leading to: 

4(n + l) Jn+im 



ak=r{-l) — 



1 + 5{m) 13 



(14) 



A particularl y important case is tha t where a(r) is proportional to the circular prolate function 
ipQfl{c,rR)^ (jSoummer et al. II2003I ). In thi s case the inte gral in Eq. (fT2]) can be computed 
using the expansion of (pQfl{c,r) derived in (jSlepian Ill964l ): 



^o,o{c, r)=Y d^jf{c)y/rF{k + 1, -k; 1; ^ 



(15) 



k=0 
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The function F( k-\-l, — fc; 1; r^) defined in Eq. (jASP reduces to a polynomial of order 2k which, 
as mentioned in ([Slepian Ill964l ) "is closely related to the Zernike polynomials" . Indeed using 
Eq. ()A7p and the results below it can be easily checked that: F{k+l, -k]l;r^) = (-l)'=i2°;.(r). 
Inserting this expansion in Eq. (|12p and integrating terms by terms leads to integrals which 
generalize Eq. ()A2p . These integrals can be c omputed for example using of integrals of the 
type Jq r''' Jm{Pr)dr ( Gradshteyn et al. 2000). This derivation will not be presented herein 
for sake of brevity. 



Finally, for more complicated compl ex amplitudes, the a^. ca n be of course computed numeri- 
cally. This problem as been addressed in (|Pawlak and Liao II2002I ) using a piecewise approximation 
of "^Aix,!/) over a lattice of squares with size A x A and centered on point {xi,yj). In this case 
the estimation of is given by: 



'{Xi,yj) 



(16) 



where 



Xi ,yi) is the integral of t he Zernike polynomial U^{p/R,(l)) over the square centered 
on (xi,yj). (jPawlak and Liao 1 12002| ) gives bound for the mean integrated squared error on the 
reconstruction of ^A{x,y) when the coefficients are given by Eq. (jlOp . This analysis is particularly 
important in our case because it quantifies the dependence of the error on the smoothness of 
^A{x,y), the sampling rate A and the geometrical error due to the circular geometry of the pupil. 



3. Coronagraph response 

3.1. Response of the coronagraph to a Zernike polynomial 

The purpose of this section is to compute the complex amplitude in D when the complex 
amplitude in A is the Zernike polynomial with radial degree n and azimuthal frequency m. In 
this case the complex amplitude ^jj{r,6) will be denoted as 'D'^{r,6). According to Eq. ([7]), the 
difficulty in the computation of T>^{r,0) lies in the evaluation of the convolution: 

H(r,^) = ((C4(r,^ + ^)m(r)) *p(r)) (r, 0) (17) 

In this expression m(r) is an "annular" mask of radius d which, with the definition adopted in 
Eq. ([3]) is defined as: 

m(r) = n (18) 

The computation of the convolution in E(r,9) is sketched in Fig. [TJ Using Eq. ([H]), H(r, ^) 
simplifies to: 

.d/Xf _ 

H(r,6i) = j J ^a(/0, (t> + 7r)p i^^/r^ + - 2rp cos{e - (/>) j pdpdcp (19) 
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The next step consists in substituting in this equation: 

• p{r) by the Fourier transform of p(r) =11 (r/R): 

, RJi(27rRr) , , 

j5 r = ^ (20) 

r 

• '^a{p^^) by U^{p/R,(l)) and consequently ^Aip,<P) by R'^U^{rR,(p) where Uj^{r,(j)) is given 
in Eq. (|A5]) . 



In order to simphfy the notations we define the new "standardized" integral H(r, 0,^) by: 



? r27T Ji ( Y^r^ + - 2rpcos{0 - (, 

{r,9,0= I I cos(mc/>)J„+i(p) \„ , ^ ' ^^^"^ (21) 

yr^ _l_ pz _ 2r pcos[t> — <p] 



10 JO 

It can be easily checked in this case that: 

j2 jn I 



E{r,e) = R'ri-l)^^ (27rRr,e,^] (22) 



2TTRd\ 



Analytical computation of H(r, ^,^) relies on the properties of the Gegenbauer polynomials 
defined in appendix [Al Substituting Eq. (jAlOj) for = 1 in Eq. (f2T]) allows indeed to separate the 
integrations with respect to p and (j): 



Eir,e,0 = 2yik + 1):^^^ [^Mp)Jn±Mdp r\osim^)ci'\cosie - <p))dcP (23) 
t^n r Jo p Jo 



,Jk+i{r) Jk+l{p)Jn+l{p) 

k=0 

• Computation of the integral on <j) is straightforward using ()A9P : 

r2Tr , ^ 



/ cos(m0)C^"'^^(cos(6' — 0))d(/) = 7rcos(m0) 5{m — k + 2q) 
Computation of the i ntegral on p relies on recursion f ormulas on indefinite integrals of products 



of Bessel functions, (jAbramowitz and Stegun Ill972l ): 



k^n, -p dp = (24) 

1'^ Jn{p) 

lo P 



rS, J f \2 1 n— 1 

/ ■-^dp = fil- MO'- MO'- UC?) (25) 

Jo P ^ 
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After computation of the integral of Eq. ([23|) . substitution of Eq. ([22]) in Eq. ^ gives the 
complex amplitude in D for a single basis function ^'^(r, ^) = U^{r/R,9): 



V-{r,9)=r-H-iy-^Rcosime) (- '^"+^(^^"^^) + , r?^,„,,(2vr^d) '^^+^(^^"^^) ] (26) 

\ fe=o / 

with fj, = R/Xf and: 

k ^ 

Vm,nMO = ik + l){^6{m-k + 2q)) [ Il^±M:lh±Mdp (27) 

ZTn Jo P 



5=0 

The corresponding complex amplitude in C for r < R can be directly computed from Eq. (j26p 
using the inverse Fourier transform of cos(m0) Jfc+i(27rr)/r obtained in Eq. ([AT]): 

C(r,0) = (-1)"^ cos(m0) (-7^- (-^) +6f;r?„,„,,(2v^Md)7^^ (^) ) (28) 

V k=0 I 

Eqs. ([26I28P give an analytical expression of the complex amplitude in C for r < R and in D 
when a single basis function is applied in A and when the size of the Lyot stop equals the size of 
the entrance pupil. In the general where the amplitude in A is given by Eqs. ([8l9p . the complex 
amplitudes in C and D become: 

^c(r,^) = «Kn)C(^,^), ^Dir,9) = J2 c^MV^^ir^O) (29) 

(m,n) (m,n) 

As mentioned in section 2.1, if the entrance pupil is not apodized a reduction of the Lyot 
stop must be considered. This is achieved replacing p{r) by p{a^^r) with a < 1. The expression 
of the complex amplitude in C is of course straightforward and for example Eq. ([28p becomes 
C'^{r,6)p{a~^r). This result allows numerical computation of the complex amplitude in D using a 
single Fourier transform. Unfortunately it is much more complicated to obtain an analytical expres- 
sion of the complex amplitude in D. The derivation presented above can be of course redeveloped 
replacing p(r) by a^p{ar) and straightforward computation shows that: 

1. Similarly to Eq. ([26p . the convolution ([17p will expand in an infinite sum of functions 
cos(m0) Jfc_i_i(27ra;ur)/r. However, the "radial contribution" to the coefficients weighting these 
functions, see Eq. ([27]) . becomes: 

Jo P 

which cannot be computed straightforwardly as in Eqs. (|24l25p . 

2. The first term in Eq. ([7]) is now replaced by the Fourier transform of Ul^ {p / R, (j))Ii{r / aR) 
which cannot be anymore calculated using Eq. (|A2p . 
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3.2. Bound for the truncation error of DJJ*(r, ( 



As we are interested in the computation of C^(r, 0) or !D™(r, ^) from the implementation of 
formula (j26p . the errors produced when the infinite sum is truncated must be studied. In order to 
reduce mathematical developments we only present herein the results for 'D^(r, 9) when the size of 
the Lyot stop equals the size of the pupil. 



We define the truncation error on T>^{r, 6): 



£^Ar(r, 9; m, n, /j,, d) = eR 



cos mt 



2^ 'nm,n,k{^'^m 



k=N+l 



(30) 



Computation of a bound on the truncation error relies on the classical upper bound for the 
Bessel functions of integer order (lAbramowitz and Stegun Ill972l ): 



|J.+i(r)|<i^^, r>0 

Substitution of this result in Eq. ([27|) gives: 

k 

Vm,n,k{0 < ik + l){Y,^im-k + 2q)) 

g=0 

k+n+2 



1 



1 



k + n + 2k\n\ \2 



k+n+2 



< 



k + l fC 
kin] 



which leads to the following bound for the truncation error: 



<?Ar(r, 9; m, n, ^, d) < 



ei?(7r^)3+"d2+" ^ k + l,. .k 



k=N+l 



(31) 

(32) 
(33) 

(34) 



The above serie is absolutely convergent for r > 0. As a consequence the expansion in Eq. (j26p 
converges uniformly for {r,9) E [0, oo) x [0, 27r). Finally, it is worthy to note that the computation 
of the infinite sum in the upper bound ([Ml) can be avoided using the equality: 



E 

fc=0 



k + l 

m 



-X 



h{2^) + ^h{2^) 



(35) 



where Iu{x) is the modified Bessel function. 



4. Simulation results 

4.1. Response of the coronagraph to the first Zernike function 

Figures [2] and [3] give the intensity in the D plane of the coronagraph when the complex 
amplitude in the A plane is one of the first six Zernike polynomials. The complex amplitudes have 
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been computed using Eq. (j26p . Each raw contains U^{r,6) and T>^{r,6) for a given couple {n,m). 
These plots have been obtained truncating the infinite summation of Eq. (j26p to the first 40 terms. 

The relevance of the truncation error bound is verified in Fig. HI This plot shows the error 
bound ()34p as a function of r for the parameters used in Figs. [2] and El The increase of the bound 
with r is simply due to the fact that the majoration of | Jfc+i(r)| given by Eq. (I3ip is only relevant 
for small values of r as long as | Jfc+i(r)| is bounded on [0, oo). It is important to note that this plot 
justifies, at least for this configuration, the validity of a truncation to A'' = 40 for the computation 
of 'Dl^{r, 9). In this case the truncation error is in fact always less than lO^^*'. 



4.2. Application to tip-tilt error analysis 



The ef fects of a tip-tilt error in Lyot cor onagraphs has been extensively studied by lLloyd and Sivaramakrishnan 



(|2005l ) and ISivaramakrishnan et al. I (|2005l ). The scope of the simulation presented here is only to 
validate the results derived in section 2 simulating the particular case where there is a misalignment 
of the star with the center of the stop. According to the previous notations the complex amplitude 
in D decomposes as Eq. ([291) . In the case of a tip-tilt error in A, the values of the coefficients a(.m^n) 
are given by Eq. (HH). 

Fig. m shows |^'£)(r, ^)| for different values of /3 > (the case /? = is given in the first row of 
Fig. [21) . The truncation in the summation ([29p has been chosen taking into account that Eq. ([14p 
implies: 

\0'{m,n)\ ~ r:^ — TTV^ I ~ I i when n ^ oo 



v/2^/?(l + (5(m)) \2nJ 

Note that according to the notations of Eq. dH), ^B{r,0) equals Eq. ([20l) shifted of -f3\f /{IttR) 
on axis x. Consequently, the star is behind the focal stop in the first two images and outside in the 
last one. 



5. Conclusion 



In this paper we have presented a theoretical formalism for the analytical study of the Lyot 
coronagraph response. The main purposes of this work are of course to assist coronagraph design but 
also to improve data processing performances for the detection and characterization of extrasolar 
planets. 



• The first application is the computation of the response of the coronagraph to a planet at a 
given position. This is achieved for example in the case of a classical Lyot coronagraph using 
Eqs. ([29ll4p . This point is essential for the derivation of an optimal decision scheme to test 
the presence of a planet at a given location. 
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This formalism can also be applied to fully characterize the statistical properties of the com- 
plex amplitude in the D plane. For a given spatial covariance in A which is fixed through the 
covariance of coefficients at, the spatial covariance in D becomes: 



coy[^> D{r,e)^D{r',e')] = y^CQy[ak,ai]Vk{r,e)Vi{r' ,e') 



k,l 



(36) 



Although detection algorithr ns based solely on th e marginal distribution of the complex ampli- 
tude can be developed as in (jFerrari et al. II2005I ). the use of an accurate model for the spatial 
correlation of the complex amplitude is ess ential in order to derive detection algorithms with 



optimal performances, as demonstrated in ([Chatelain et al. II2006I ) 



The author thanks the anonymous referee who helped improve the paper. The author is also 
grateful to Claude Aime and Remi Soummer for helpful discussions and insightful comments. 



A. Appendix 



This section presents some facts about Fourier transform in polar coordinates, Zernike and 
Gegenbauer polynomials. 

Among the various available possibilities to define an orthogonal set of funct ions on the unit 
radius disk a central position is hold by the Zernike polynomials, see for example (jMahaian Ill994l ) 
and included references. They are defined for n > m by: 



R^{r)cos{me)Il{r), U'^'ir^e) = R'^{r)sm{me)U{r) 



(Al) 



when n et m share the same parity. The R^{r) ar e the radial polynomia ls. Different normalizations 
exist for R^{r), we retain herein the definition of (jBorn and Wolf Ill99ll ): R^{1) = 1. Among many 
properties verified by these polynomials, we focus on: 



(A2) 



see (iBorn and Wolf I Il99ll . appendix VII) for the proof. This equality allows straightforward 
computation of the Fourier transform of the Zernike polynomials. In fact recall first that when 
/(r, 9) = g{r) cos(m0), m G Z, a simple change of variables in the Fourier transform integral leads 
to: 

(A3) 



(A4) 



/(p, (/)) = 27r(-j)™'cos(m0) / rg{r)Jmi2TTrp)dr 

Jo 

An analog result for the inverse Fourier transform of f{p, 0) = h{p) cos(m</>) is: 



f{r,e) = 27rj™cos(me) / ph{p)J^{2Tirp)dp 

Jo 
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Applying the result of Eq. ()A3P with Eq. (|A2p immediately gives: 



P 



= j™(-l) 2 sm{m(l)) ^ 



(A5) 
(A6) 



The previous equation gives the inverse Fourier transform of cos{m(j)) -^"+1^^^^) when n > m > 
and n et m share the same parity. In the general case where n > and m > this inverse Fourier 
transform, denoted as /(r, 6) must be computed independently. If we subs itute h(r) by Jn+^ (2Trp)/p 



in Eg . ()A4p the resulting integral is a Weber-Schafheitlin type integral (jAbramowitz and Stegun 
I972I ). This results in /(r, 6) = j*" cos(m6l)7^™(r) where: 



ifr<l, n^{r) 



r(^^ + i) 



r(m + i)r(M + i) 



F 



n + m 



+ 1, 



m — n 



-; m + 1, r" 



(A7) 



(A8) 



F{a,b; c; z) is the Gauss hypergeometric function, see (jGradshtevn et al. Il200d ): 

, , ab a(a + l)b(b + l) o 

It is interesting to note from Eqs. ()A7p and ()A8P that if 6 = (m — n)/2 € Z~ the sum in Eq. (jASP 
reduces to a polynom in z of order — (m — n)/2. Consequently lV^{r) reduces to a polynom with 
degree n which of course coincides up to (— 1)(™~")/2 with R^{r) for r < 1. For this reason TZ'^{r) 
can be considered as a natural ge neralization o f the Zernike polyno mials. Note that, contrarily 
to the generalization proposed in (jMyrick 1119661 ) or (jWiinsche I l2005l ) . this generalization is not a 
polynomial. 

We i iow briefly give the p rinci pal results related to the Gege nbauer polynomials. See for 
example ( Andrews et al. Ill999l ) or ( Abramowitz and Stegun I Il972l ) for detailed properties. The 
Gegenbauer (or ultraspherical) polynomials, noted as i 1-^ cj!^\t) are defined as the coefficients of 
the power series expansion of r (1 — 2rt + r^)""^: 

^ 00 



(1 -2rt + r 



2\u 



k=0 



For example cj^\t) gives the Chebyshev polynomial of the second kind Uk{t): 



Cj^\cosm = Y,cos{{k-2q)^l;) 

q=0 



(A9) 



Among the numerous beautiful properties of the Gegenbauer polynomials, we focus on the 
expansion: 

(AlO) 



= 2^r{u) f^{k + .)^%^%^cf (cos(7)) 



k=0 



where w = y^r^ + — 2r/9 003(7). 
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This preprint was prepared with the AAS lAT^^X macros v5.2. 
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Fig. 1. — Computation of the convolution between ^Ai^, Q + vr)m(r) and p(r). 
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Fig. 2. — Complex amplitude in A and squared root of the amplitude in D, i.e. \'D^{r,9)\. The 
parameters used in the simulation are: Xf = 1, R = 1, d = 3, e = 1 (Lyot coronagraph) . 




Fig. 3. — Complex amplitude in A and squared root of the amplitude in D, i.e. \T>J^{r,9)\. The 
parameters used in the simulation are: Xf = 1, R = 1, d = 3, e = 1 (Lyot coronagraph) . 
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Fig. 4. — Bounds on the truncation error as a function of r. The parameters are the same as the 
parameters used for Figs. [2]and[3l For each value of N, the bound is plot for the first 6 Zernike 
polynomials. 
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Fig. 5. — |^z)(r, 0)1 for different values of /3. The parameters are the same as the parameters used 
for Figs. El and El 



